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ABSTRACT 


A thermal analysis of Hg/Cd/Te solidification in a Bridgman cell is made using 
Continuum’s VAST code. The energy equation is solved in an axisym metric, quasi-steady 
domain for both the molten and solid alloy regions. Alloy composition is calculated by a 
simplified one-dimensional model to estimate its effect on melt thermal conductivity 
and, consequently, on the temperature field within the cell. Solidification is assumed to 
occur at a fixed temperature of 979 K. Simiplified boundary conditions are included to 
model both the radiant and conductive heat exchange between the furnace walls and the 
alloy. Calculations are performed to show how the steady-state isotherms are affected 
byj a) the hot and cold furnace temperatures, b) boundary condition parameters and c) 
the growth rate which affects the calculated alloy’s composition. 

The Advanced Automatic Directional Solidification Furnace (AADSF), developed by 
NASA, is also thermally analyzed using the CINDA code. The objective of this analysis is 
to determine the performance and the overall power requirements for different furnace 
designs. 
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INTRODUCTION 

This contract (NAS8-35331) started on June 17, 1983 with an intended period of 
performance of 12 months. The satisfactory progress of the research resulted in a series 
of modifications primarly intended to increase the period of performance and streamline 
the scope of work as necessitated by technical developments. The modified program can 
be simply divided into two main phases. The results and conclusions from both phases 
have already been delivered to NASA. 

The first phase was the thermal analysis of the growing crystal inside the Bridgman 
cell. This phase was completed in March, 1985'. An interim report was submitted in 
August 1984. A copy of this interim report is attached as Appendix A. The thermal 
analysis of the growing crystal is now being continued and expanded under a separate 
contract, NAS8-36483. 

The second phase started in March, 1985 and was concerned with thermal analysis of the 
solidification furnace, its performance, and overall power requirements. Different 
furnace designs and operating conditions were supplied by NASA and analyzed using the 
CINDA code. The results of these calculations were continuously delivered to NASA so 
that the furnace design process could continue uninterrupted. No father documentation 
of these results is required since they have already been delivered to the satisfaction of 
NASA's contracting officer representatives. 
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TECHNICAL DISCUSSION 

Mercury Cadmium Telluride crystals are essential for manufacturing infrared detection 
devices. The purity and composition uniformity of the crystals are most important, and 
unfortunately, difficult to achieve. The objective of this work is two fold: 1) to 

calculate the temperature distribution within the sample itself during solidification under 
different operating conditions and, 2) to thermally analyze the directional solidification 
furnace and determine it's preformance and power requirements. 

The first phase of the investigation specifically addresses the thermal phenomenon which 
dominates the process of crystal growth. Numerical axisymmetric heat conduction 
within the Hg/Cd/Te sample is performed. The sample is regarded as consisting of a 
solid zone and a melt zone. The thermal conductivity in the melt zone is temperature 
and composition dependent. The interface between the two zones is taken as the 706 C 
isotherm. The composition in the solid zone is considered uniform and equal to the initial 
bulk composition. The heat exchange between the furnace and the sample through the 
air gap and the glass ampule is a complex phenomenon of radiation and conduction. This 
heat exchange is approximately calculated by a simple radiation model and a heat 
balance on the glass wall elements. The resulting energy flux at the sample's surface is 
then used as a boundary condition for solving the heat conduction equation within the 
sample's domain. The solution provides the thermal field within the sample and, hence, 
the shape of the interface. Parametric cases are calculted to study the effect of the 
furnace operating conditions on the crystal's temperature field and interface shape. 
Complete details of the calculations, results, and conclusions are presented in Appendix 
A as a reproduction of the previously submitted interim report. This work is now being 
expanded and continued under a new and separate contract (NAS8-36483). In this new 
effort, an accurate radiation model is used to calculate the energy exchange between the 
furnace and the sample's surface. Also, the composition and conductivity of the solid 
phase are more realistically determined as functions of the interface shape and the 
temperature distribution. The interface itself is no longer an isotherm, but rather is 
obtained from the empirical information of the Hg/Cd/Te phase diagram. The slow 
circulating motion in the liquid zone is also considered. 

The second phase of the investigation specifically addresses the performance of the 
directional solidification furnace. The CINDA code is used to calculate the thermal field 
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within the components of the furnace and its overall power requirements. This is an 
ongoing study in support of the furance design effort carried out by NASA. The furnace 
geometry and construction materials are changed and updated in search of an optimum 
design. Continuum supplied corresponding predictions to guide the designers. All results 
and calculations were continuously delivered to the satisfaction of NASA's contracting 
officer representatives. Thus, no further documentation is required. 
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ABSTRACT 


A thermal analysis of Hg/Cd/Te solidification in a Bridgman cell is made using 
Continuum’s VAST code. The energy equation is solved in an axisym metric, quasi-steady 
domain for both the molten and solid alloy regions. Alloy composition is calculated by a 
simplified one-dimensional model to estimate its effect on melt thermal conductivity 
and, consequently, on the temperature field within the cell. Solidification is assumed to 
occur at a fixed temperature of 979 °K. Simiplified boundary conditions are included to 
model both the radiant and conductive heat exchange between the furnace walls and the 
alloy. Calculations are performed to show how the steady-state isotherms are affected 
by: a) the hot and cold furnace temperatures, b) boundary condition parameters and c) 
the growth rate which affects the calculated alloy's composition. 
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INTRODUCTION 


Bridgman-type crystal growth techniques are attractive methods for producing homo- 
geneous, high-quality infrared detector and junction device materials. However, crystal 
imperfections and interface shapes still must be controlled through modification of the 
temperature and concentration gradients created during solidification. The objective of 
this investigation was to study the temperature fields generated by various cell and heat- 
pipe configurations and operating conditions. Continuum's numerical model of the 
temperature, species concentrations, and velocity fields (1) was used to describe the 
thermal characteristics of Bridgman cell operation. Detailed analyses of the other 
transport processes will be addressed in subsequent phases of this research program. 

Alloy in the Bridgman cell is melted by an enclosing furnace (a high temperature heat- 
pipe); the alloy is solidified by slowly replacing the heated walls with cooler walls. 
Radiative exchange between the furnace walls and the alloy is a primary means of heat 
transfer, as is the conduction of heat along the glass ampule wall. Tailoring of cell 
operation by control of the heat transfer process is described qualitatively in (2). This 
study will address cell design by quantitative evaluation of the thermal features of the 
system. 

Experimental studies by Lehoczky and Szofran of Cd/Hg/Te alloys in Bridgman cells (3) 
and of the phase equilibrium behavior of these alloys (4) indicate that an HgTe enriched 
layer should accumulate on the melt side of the interface, hence both thermal and solutal 
density gradients would tend to stablize the fluid and prevent buoyant convection. 
Furthermore, local measurements of crystal composition showed significant radial 
variation (3). Measured axial composition variations were interpreted in terms of phase 
equilibrium relationships and axial diffusion, but the radial variations were attributed to 
other factors. The implication is that radial temperature gradients must affect the 
shape of the melt/solid interface and, thus, create thermal and/or solutal currents. This 
study predicted melt/solid interface and isotherm shapes for a variety of wall heating 
conditions. Velocity and species distribution effects were not considered in this work, 
although momentum and species conservation equations are already included in the 
Bridgman cell 'model. The very complex heat transfer boundary conditions for the cell 
must be accurately modeled before it is meaningful to attempt predictions of the 
convective flow caused by curved interface shapes. 
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The geometry shown in Figure 1 is complex and does not lend it'self to parametric 
analyses. Rather cell and furnace dimensions reported in reference (3) will be used 
exclusively and parametrics were reserved for surface radiative properties and glass 
transmissivity, which are the least well known of the important thermal properties of the 
system. Effects controlled by system operating conditions such as the hot and cold heat 
pipe temperatures and growth rates were investigated. 

The analyses performed and the predicted results are presented in the remainder of this 
report. Critical design conditions were identified for further study in the later phases of 
this research program. 


TECHNICAL DISCUSSION 

A thermal analysis of the Bridgman cell was made for the axisymmetric, quasi-steady 
system shown in Figure 1. The energy equation was solved for both the molten and solid 
alloy regions. Solidification was assumed to occur at a fixed temperature. One-dimen- 
sional species solutions were used to estimate the effect of composition on melt thermal 
conductivity, otherwise the energy equation solution from Continuum's Bridgman cell 
model was used. Boundary conditions along the glass wall are complex due to (a) the 
necessity of treating the glass with a lumped parameter model (or as a separate 
computational region) to determine the temperature distribution within the glass with 
sufficient accuracy to evaluate the heat conduction to the alloy and (b) the integral 
nature of the radiation flux in the annular space between the furnace and the alloy 
requiring a major assumption (or a three-dimensional view factor evaluation and an axial 
integration of the heat flux). The simplified flux equation was used in this study to 
expedite generating a wider range of parametric results for an initial thermal analysis. 
Details of the calculations and a description of the alloy properties used are given in 
subsequent paragraphs; followed by a presentation of parametric cases which were 
studied. Two types of analyses were made: in the first, simple heat conduction within 
the cell was predicted for specified wall temperatures to compare to solutions which 
have been reported by Naumann and Lehoczky (5) and, in the second, various wall heating 
conditions wefe studied. 
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Continuum's VAST code was used as the basic Bridgman cell model (1). Options of 
including radiation and conduction to the cell walls, as well as conduction and transmis- 
sion through the fused silica walls of the ampule are available. The top and bottom ends 
of the ampule are too far removed from the interface to effect the temperature field in 
the vicinity of the interface, therefore temperature boundary conditions are appropriate 
at these locations. The VAST code is a numerical solution to the species, momentum, 
and energy equations. Die analog utilizes an explicit time integration to produce a 
steady-state solution for a nodal network which is obtained from a finite-element formu- 
lation. For this study, only the energy equation solution was used. The Bridgman cell 
model contains a detailed specification of boundary conditions to describe heating 
through the glass wall and from the furnace walls; these boundary conditions directly 
couple the environmental conditions of the alloy in the cell to the node temperatures 
along the outermost edge of the alloy. The only time effect considered was to specify 
the molten alloy thermal conductivity as a function of growth rate, therefore the only 
significant thermal effect omitted from this model is the heat of fusion which is a direct 
function of growth rate. A more accurate analysis cannot be obtained without also 
solving the species equations. 


PHYSICAL PROPERTIES OF THE Hg/Cd/Te SYSTEM 

Best available physical property data for the Hg/Cd/Te system are summarized in (1); 
however, to expedite an initial thermal analysis study, some simplification and empirical 
curve fits of these data were used. The following alloy properties were used for this 
study. The density of Hgj_ x Cc^ Te alloy was assumed constant in all calculations and 
equal to 7.54 gm/cm 3 . This corresponds to a liquid phase mole fraction, x, of 0.2 which 
is the initial composition. Density variations are neglected since convective motions 
inside the ampule are not considered. 

The local alloy composition, x, is calculated at each node in the liquid phase using the 
following one-dimensional species solution (6): 

x = x„ jl - exp( -RZ/D) | (for Z > 0) 
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where S is a segregation coefficient equal to 3.9; R is the growth rate; D is an effective 
diffusion coefficient equal to 5.5 E-5 em^/sec; and Z is the vertical distance above the 
melting isotherm which is easily determined from the temperature distribution. 
According to the above equation, x is equal tox 0 /S at Z = 0 and changes 
exponentially in Z to approach x 0 as Z becomes large. In the solid phase (Z < 0), the 
composition is assumed constant and equal to x o . 

The thermal conductivity, K, is calculated at each node in the liquid phase using the 
following equations (5): 

a = B ( ln(T) ) - A 

K = a p C p 


where a is the thermal diffusivity; T is temperature in C; p is the constant density (7.54 
gm/cm^); Cp is the heat capacity (1.8428 E 6 em^/sec^ °K). The empirical coefficients A 
and B are given in Table 1 as functions of composition, x. In the solid phase, the thermal 
conductivity is assumed constant and equal to 6.3947 E 4 gm cm /sec 0 K. The heat of 
fusion was neglected in this analysis. 


Table 1. Values of the Empirical Coefficients A and B. Data from Ref. (6) 


X 

A 

B 

0 . 

67.401 

10.4655 

0.052 

46.767 

7.25493 

0.107 

46.024 

7.09724 

0.205 

42.041 

6.44051 

0.301 

37.423 

5.69691 
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BOUNDARY CONDITIONS FOR A BRIDGMAN CELL 


Thermal energy is exchanged between the alloy and the furnace, through the glass wall of 
the ampule, by radiation and conduction. A part of the radiated energy is directly trans- 
mitted to the alloy while the glass wall absorbs the remaining part and reflects a negli- 
gible amount. The radiant energy absorbed by the glass, along with the energy conducted 
through the air gap, is conducted through the two-dimensional domain of the glass wall. 
A part of this conducted energy reaches the alloy, while the remaining part flows in the 
longitudinal direction. The net result is a certain amount of energy reaching or leaving 
the outer surface of the alloy per unit area per unit time. An accurate evaluation of this 
energy flux is vital for predicting a correct temperature distribution within the alloy. 
Once calculated, this energy flux, Q, can easily be used to impose the proper boundary 
condition on the solution of the energy equation. 


In order to calculate the boundary heat flux, Q, a "lumped” heat balance is performed on 
the portion of the glass wall associated with each boundary node (see Figure 2). This 
element of glass is assumed to have a uniform temperature T a which is an average of the 
node temperature, T n , and the temperature of the outer surface of the ampule at a point 
facing the node, Tg. T n is obtained from the solution of the energy equation at the 
previous time iteration and Tg is calculated by equating the heat conducted to and 
through the air /glass interface: 


<K a /5,) <V T g ) = <K g /S g ) <T g - T n ) 


where K. and K„ are the conductivities of air and class: 6 and $ are widths of the air 

“ O fl 2 

gap and glass wall; and Tf is the furnace temperature at a point facing the node. In a 
similar manner, T ga and T ab are calculated to represent average temperatures in the 
glass portions above and below the portion under consideration. 

The set of temperatures Tf, Tg, T ft , T n , T aa and T &tj is used to calculate the following 
heat fluxes: 




a 


fir / o \ 

'V°za' 


/ rr 

V 1 


aa 


\ 

/ 


«h = V 4 cb> (T .b- V 
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S = (J VV ( V v 


Q = <JK /S > <T - T > 
n g g n a 


Q r = H r <T f '- T,*) 


Q (absorbed) = Q (1 - t ) 
r r 

Q p ( transmi tted) = Q r x 


where t is the transmissivity of the glass, and H p is a radiation coefficient the value of 
which is discussed later. In each of the above formulae, Kg is the glass conductivity at a 
temperature equal to the average of the pair of temperatures used to calculate the 
flux. Also, t is the transmissivity of the glass at a temperature equal to (Tf + T fi )/2. 
Figure 3 shows the temperature dependency of K a , Kg and t . 

The energy fluxes cause the temperature of the glass element to change from T„ to T 

w u U 

The new T a is calculated from the heat balance on the glass ’lump’': 


(V g p g C g / At ) (T; - T a ) = (Q a + Q b ) Aj + (Qg+ Q p ( absorbed) )A^ + Q n A 2 
where V_ is the volume of the glass element; p is the density of the glass 

® a D 

(2.203gm/cm°); Cg is the specific heat of glass; and A t is the time step size. The areas 
Aj, Ag and Ag are shown in Figure 2. Finally, the flux Q which represents the energy 
exchange between the furnace and the outer surface of the alloy in the vicinity of the 
node is calculated as: 

Q = (W tr /6 tr ) <T’ - T ) + Q (transmitted) 

K g ti n r 


At each boundary node, the energy per unit volume p (obtained by solving the 
energy equation) is adjusted to accommodate the flux Q as follows: 

pE new * pE old + « Al 4t/v 
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Where v is the volume of the alloy element associated with the node. This energy adjust- 
ment is performed at the end of each time step. 

In calculating the radiative heat flux, Q r , a radiation coefficient, H r is used. To esti- 
mate H r , the following assumptions are made: 

1. The area A 2 (see Figure 2) is a gray surface with a uniform temperature of T fl . 

2. The area A 2 "sees" a furnace area, Af, which is also a gray surface at a uni- 
form temperature of Tf. This is a reasonable approximation since A 2 sees 
primarily that part of the furnace radially facing it. Tf is taken to be the hot 
furnace temperature for all nodes above the barrier, and the cold furnace 
temperature for all nodes below the barrier. For the nodes facing the barrier, 
Tj varies linearly between the hot and cold temperatures. 

/ 

3. Both A 2 and Af satisfy Kirchhoff's Law, i.e. the absorptivity, a, is equal to the 
emissivity, e - 

Based on these assumptions, the idealized gray-body heat balance (7) is used: 

Q r = o <T ( '- T a ‘) /|(Aj/A f > (p { / Ef ) + 1/F + (P al /c al >) 

where o is the Stefan-Boltzmann constant (5.6696 x lCT^gm/sec^ °K^); p ^ and p fll 
are reflectivities; e ^ and j are emissivities of the furnace and the alloy’s surface 
respectively. F is a view factor between A 2 and Af. Since A 2 "sees" Af and nothing 
else, F is equal to 1. Furthermore, the furnace area Af, which is much larger than the 
elemental area A 2 , has a large emissivity and a small reflectivity. Thus, the term 
(A 2 /Af) ( Pf / Cf ) is neglected and an approximation of H r is given by: 


Using c , = 0.2 and p . = 0.8, the estimated value of H becomes 1.1 x 10”® 
(gm/sec J °K 4 ). The value of H f obviously controls the efficiency of the radiative 
exchange between the furnace and the alloy. Using H f as a parameter should demon- 
strate the impact of the radiation mode of heat transfer on the calculated temperature 
field within the alloy. 

The radiative energy exchange between the furnace and the alloy can be modeled more 
precisely by an expression for the net energy radiated from a differential area on the 
furnace to a differential area on the alloy's surface per unit time. Intergrating this 
expression over Aj and over the furnace area "seen" by it gives an accurate value of the 
radiant energy reaching A 2 per unit of time. The determination of whether or not this 
method appreciably improves the accuracy of the radiation model is the subject of a 
subsequent investigation by Continuum. 


PARAMETRIC STUDIES 


The model described above was solved, using Continuum's work station facilities, for a 
wide range of input parameters. The calculations are divided into four groups of runs. 
The first group is intended to verify the performance of the model by requiring it to 
duplicate an appropriate analytical solution. Each of the last three groups is designed to 
investigate the effect of a certain input parameter on the alloy's temperature distribu- 
tion at the steady state. 

The first group consists of two runs. In Run #1, the case of infinite solid cylinder, with 
constant conductivity and a given surface temperature, is simulated. Therefore, neither 
the composition calculations nor the boundary conditions described above are invoked, 
but the energy equation is numericaly solved in the axisymmetric domain with a specified 
set of temperatures at the surface nodes. Run #2 differs from Run #1 in that the 
conductivity of the material (tops from 7K to K at all locations where the temperature 
is less that 700 °C . The specific value of K has no effect on the results. Figure 4 shows 
a comparison between the analytical solution reported in reference (5) and Continuum's 
solution. ' 
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In the second group of runs (Runs #3, 4, 5), the effect of the furnace temperature 
distribution is investigated. Composition and boundary condition calculations, as 
described in the preceding sections, are invoked. The melting point temperature is 979 

0 

K and the growth rate is 6 x 10 cm/sec. The rest of the parameters are as mentioned 
earlier. The hot and cold furnace temperatures, and the steady state temperature 
contours within the alloy are given in Figure 5 . It is clear that the furnace temperatures 
have a great effect on the location and curvature of the alloy's isotherms. 

The third group of runs (Runs #6, 3, 7) shows how the alloy's steady state isotherms are 

affected by the rate at which heat is exchanged between the alloy and the furnace. This 

rate is controlled by the magnitude of the radiation coefficient, H r , and the width of the 

air gap, . A smaller H r means less radiative exchange while a smaller 6 fi causes 

larger conductive exchange to occur. The values of H r and for the three runs are 

given in Figure 6 and indicate that energy exchange is most efficient in Run #6 and least 

efficient in Run #7. Run #3 uses the real value of 6 and the value of H_ estimated in 

a r 

the preceding section. Concerning the "efficiency” of heat exchange at the boundary, 
Run #3 is intermediate between Runs #6 and 7. For all of the three runs, the hot and 
cold furnace temperatures are 1099 °K and 788 °K respectively. The growth rate is 6 x 
10 -6 cm /sec and the rest of the input parameters are as mentioned earlier. The steady 
state temperature distributions are shown in Figure 6 and compared to the constant 
temperature distribution of the furnace. As expected, the wall temperatures of Run #6 
are closer to the furnace temperatures than those of Runs #3 and 7. The temperatures 
at the cylinder's axis show the same behavior but to a lesser extent. 

The last group of runs (Runs #4, 8, 9) shows the effect of the growth rate on the alloy's 
temperature balance. In all the three runs, the hot furnace temperature is 1273 °K and 
the cold temperature is 673 °K . The radiation coefficient is 1079 x 10 - ® gm/sec® °K* 
and the width of the air gap is 11.5 mm. The growth rates of the three runs are given in 
Figure 7 along with the corresponding temperature contours at the steady state. 
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DISCUSSION OF RESULTS 


The comparisons shown in Figure 4 demonstrate the accuracy of the numerical model 
used to simulate the Bricjgman cell. The modest grid density used for these calculations 
causes no significant error in the results. Notice that no heat flux is predicted across the 
cell centerline, although no explicit specification of zero gradients is imposed on the 
centerline. 

O 

The cell model used assumes that solidification occurs on the 979 K isotherm. The hot 
and cold temperatures which are chosen for furnace operation move the melting 
isotherm, and, as shown in Figure 5, the shape of the isotherm is determined by where it 
is positioned in the celL The heat transfer boundary conditions obviously control the 
isotherm location. Even though the individual boundary condition specification can be 
made more accurately, it is apparent that the interface shape can be controlled by 
furnace operating conditions. For the radiation boundary conditions and system physical 
properties used, Figure 6 shows that the interface is only slightly affected by the radi- 
ation and conduction levels. A more precise radiation analysis in the vicinity of the heat 
barrier region and a more accurate radiation/conduction analysis in the glass wall may 
affect the interface shape; such studies are in progress, but no results are available at 
this time. 

The relationship between interface shape and growth rate is shown in Figure 7. This 
effect is small, even though isotherms away from the interface are markedly affected by 
the growth rate. Exact heat transfer boundary conditions in the barrier zone region will 
have some effect on the growth rate dependence of the interface, but the major effect is 
the neglected heat of fusion which is directly proportional to the growth rate. The 
experimental studies reported in (3) show flat interfaces for Hg/C d /T e systems for very 
slow growth rates and curved interfaces for all of the larger growth rates conditions. 
The effect of composition alone does not appear to explain radial variations of interface 
shape because conductivity changes only slowly from the solid value at near interface 

locations in the liquid. Further analyses of both composition and thermal effects which 

* / 

occur on solidification are the subject of ongoing research. 
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CONCLUSIONS AND RECOMMENDATIONS 


Thermal modeling studies of the Bridgman cell allow us to make the following conclu- 
sions. 

1. An accurate thermal model of the Bridgman cell has been developed. 

2. Less complete thermal models of the Bridgman cell, even if they include con- 
vective flows, cannot provide accurate simulations. 

3. More detailed analyses of radiation effects in the vicinity of the thermal bar- 
rier and of the temperature distribution within the glass wall, along with better 
furnace wall and alloy surface radiation data, would improve the accuracy of 
the thermal model. 

4. A better description of the concentration gradients and heat of fusion releases 
which accompany solidification would improve the thermal analysis. 

5. The interface shape control methods reported in (2) can be quantitatively 
evaluated with Continuum’s Bridgman cell model. 

It is recommended that: 

1. The radiation analysis at the barrier and the conduction analysis within the 
glass wall be extended to include details with regard to view angles and conduc- 
tive node points in the glass. 

2. The species conservation equation with attendant interface species and energy 
boundary conditions be included in the thermal model. 

' f 

3. A wider range o? parametric studies be made. 

4. The thermal and solutal convection effects be included in the thermal analysis 
of the Bridgman cell. 
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